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Abstract 



We study the decay B — > Xg^ within the framework of the minimal 331 model, 
taking into account both new experimental and theoretical developments that allow 
us to update and improve on an existing ten year old analysis. In contrast to several 
other flavor changing observables that are modified already at tree level from a new 
Z' gauge boson, we have only one loop contributions in this case. Nevertheless, 
these are interesting, as they may be enhanced and can shed light on the charged 
gauge boson and Higgs sector of the model. Numerically, we find that the Higgs 
sector, which is well approximated by a 2 Higgs doublet model (2HDM), dominates, 
since the gauge contributions are already very strongly constrained. With respect 
to -B ^ Xg^, the signal of the minimal 331 model is therefore nearly identical to 
the 2HDM one, which allows us to obtain a lower bound on the charged Higgs 
mass. Further, we observe, in analogy to the 2HDM model, that the branching 
fraction can be rather strongly increased for small values of tan f3. Also, we find 
that B Xg^ has no impact on the bounds obtained on rare K and B decays in 
an earlier analysis. 



1 Introduction 



Presently, all phenomena observed in nature (with the exception of gravity) are described 
within the standard model (SM) of particle physics. This model seems to work beautifully 
up to scales of at least O(TeV). On the other hand, successful as the SM may be, the 
general hope (and belief) of particle physicists is that it should be a part of a more 
fundamental theory. This belief is based on several theoretical shortcomings of the 
model, the most important of which are the instability of the Higgs mass as well as the 
general particle content of the model, such as the fact that each fermion type appears in 
three generations. This last point is addressed in the context of 331 models [1,2], where 
the requirement of anomaly cancellation combined with the asymptotic freedom of QCD 
force the number of generations to be exactly three. To achieve this, the electroweak 
SU{2)l of the SM is extended to an SU{3)l, where the third generation is treated 
differently, i.e. is transformed as an anti-triplet. 

This set-up leads to the existence of several new particles, in particular of new gauge 
bosons, such as a Z' that transmits flavor changes at tree level due to the different treat- 
ment of the third generation. As a consequence, this gives rise to tree level contributions 
for several observables that proceed at loop level only in the SM and which have been 
extensively studied in the literature [3-7]. In addition to these observables, it is also 
interesting to investigate the inclusive decay B — s> ^^7, where Xg denotes a sum over 
all final states containing a strange quark. While tree level contributions are absent 
here, those at one loop may be interesting, and have first been studied for the minimal 
model (to which we will also restrict ourselves) in [8], after similar, analogous effects 
to e'/e have been investigated in [9]. There are several reasons that make this analysis 
worthwhile, in spite of the fact that a first reasoning would suggest to focus on tree level 
Z' exchange only. These are: 

• The mechanism that causes FCNCs at tree level (i.e. the different treatment of the 
third generation) also enhances the one loop contribution rather strongly. This is 
due to the breakdown of the GIM mechanism, which acts very effectively in the 
SM. In any case, it was found in [8] that the new contributions could, in principle, 
be almost comparable in size to the SM. Since the one loop effects are also governed 
by other particles than the tree level ones, i.e. charged gauge bosons and Higgs 
fields, one may also hope to shed some light on an entirely different sector of the 
model. 

• In the ten years since the analysis of [8], the data on B ^ Xs7 have improved 
considerably [10-12], leading to a very precise experimental number. This can 
be used to place constraints on the parameter space of the model which can be 
compared with those coming from other FCNC processes, obtained for example 
in [6]. 

In view of this, we update the B — > Xg'-f analysis of [8] using new data, retaining the 
general scheme of that analysis. This concerns in particular the treatment of the QCD 
corrections via the renormalization group equations, and the treatment of the Higgs 
sector. Here, we study only an effective 2 Higgs doublet model (2HDM) type II, which is 
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a good approximation in the flavor sector [3,8]. This analysis extends and concludes our 
study of FCNC processes begun in [6]. Our paper is then organized as follows: First, in 
Section [2], we introduce very briefly the minimal 331 model, focusing on those parts of 
the model that are important for the penguin diagrams contributing. Next, we give the 
most general background on B ^ Xs'j in Section [3l where we also list and explain the 
different contributions to the decay and discuss the cancellation of divergences. Section 
m then contains all our numerics, where we compare these new constraints with those 
from the measurement of the mass difference AM,,. These are the only relevant ones 
in this context, since the different sectors of flavor transitions sd, hd and hs decouple. 
Finally, we conclude in Section [5l 

2 The Minimal 331 Model 

The minimal 331 model has been discussed extensively in [3, 13, 14] and in short in 
[6], from where we take the conventions. Many variations of this minimal model have 
been developed, such as ones with right handed neutrinos, supersymmetric versions, and 
others [15-28]. Generally, in 331 models the electroweak SU{2)l of the SM is extended 
to a SU{3)l, which is broken down in two steps: 

SU{3)c X SU{3)l X U{l)x ^ SU{3)c x SU{2)l x U{1)y SU{3)c x f/(l)e™ (1) 

Evidently, this breaking process involves three Higgs doublets, one of which develops 
a vacuum expectation value (VEV) at a high scale, while the other two VEVs are of 
the order of the weak scale. This leads to a much richer scalar sector than in the SM, 
which has been extensively studied for the minimal model in [29,30], as well as a more 
complicated Yukawa structure [3,4,13,14,19]. Physically, after absorbing the appropriate 
number of Goldstone bosons, one is left with one light neutral Higgs scalar, 7 heavier 
neutral Higgs particles as well as 4 singly and 3 doubly charged Higgs particles, all of 
which are heavy. On the other hand, the quark doublets of the SM are extended to 
triplets by adding an additional heavy quark. All of the quarks, as well as their right- 
handed counterparts, have three quantum numbers called X, T3 and Tg, corresponding 
to the diagonal generators of SU{3). 

From these, the electric charge can be obtained by 

Q = n + f3Ts + X; (3 = V3 (2) 

in our normalization of the hypercharge X. To ensure anomaly cancellation, the third 
generation couples differently than the first two. This has the important consequence of 
generating flavor changing neutral currents at tree level. Another curiosity of the minimal 
331 model is the Landau Pole that can arise at rather low energies. It becomes apparent 
when one expresses the ratio of SU (3) and SU (2) couplings through the Weinberg angle 
as 

gj^ ^ 6 sin^ Ow .^x 
1-4 sin^^iy ' 

Clearly, the theory is ill defined if sin^ 9w grows to be 1/ 4, which puts an upper limit 
on the scale of the symmetry breaking as well as on the heavy gauge boson masses. 
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Analyzing this carefully [31], gives an upper bound of several TeV for the Z' mass, so that 
in the numerical analysis below we will follow [6] and use values of up to M^/ = 5 TeV. 
In choosing /? = a/3 in ([2]), and assigning the appropriate hypercharge structure, we have 
distinguished the model to be the minimal 331 model. 

Turning now to the gauge boson content, one finds that the breaking process sketched 
above produces, in addition to the SM gauge bosons, an additional neutral Z' boson as 
well as singly charged and doubly charged bosons. The photon is massless, 
as required, since a UiX) remains unbroken while the W and Z masses are at the weak 
scale, as required. Finally, all additional gauge bosons obtain their masses from the large 
VEV and are therefore heavier. Assuming that this heavy VEV is much larger than 
the others, an interesting relation between the heavy gauge boson masses can be found, 
which allows to express the Y^ mass in terms of the Z' mass or vice versa: 



_ 3 (l-4sin e^) 

4 (cos2 6',4/) ^ 

while My± = My±±. We will be more explicitely concerned with the respective fermion 
couplings below, but let us already here state that in the quark sector the new charged 
gauge bosons always transmit a coupling from a SM light quark to an additional heav 
one, and therefore do not modify the amplitudes of low energy observables at tree levet. 
The neutral Z', on the other hand, also has couplings to light quarks only, and does give 
this kind of contributions. In the decay B ^sl^ however, both charged and neutral 
gauge bosons are equally important, so that information on the charged bosons may be 
obtained. 

After these more general remarks, let us now show explicitely the respective fermion 
couplings. Concerning the neutral currents, the Lagrangian for the above-mentioned 
FCNCs at tree level is given by 

Lpcnc = ^j^^=^'^1^1lUI[ \uj^u^~d^^^jyl[ ]vLd]Z'r (5) 

The explicit couplings to fermions and a possible parameterization for the new mixing 
matrices Ul and Vl, that diagonalize the up and down-type Yukawa couplings, respec- 
tively, have been given in [6]. These obey 

UIVl = VcKM, (6) 

and the notation of ([5]) and IQ is to be understood as in [6], so that the tilde distinguishes 
between the SM CKM matrix Vckm and the mixing matrix for the down-type quarks. 
Additionally, the charged current vertices in this basis are then 

= uriLUlVd = uY-fLVcKMd 





^This is no longer true in the lepton sector, where for example a new tree level diagram to muon 
decay modifies the Fermi coupling constant [6]. 
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Figure 1: New magnetic penguin diagrams contributing to the decay B Xs'y. V denotes the heavy 
charged gauge bosons , . We do not show self-energies and the diagrams involving Goldstone 
bosons. 



= d^^LV^ \0l \ D + T^^'^L (0 1) Ulu 




J^^^ = uYilUI 1 Ld - T^^L (0 Ol)Vd. (7) 



From these equations, it becomes obvious that there can be new penguin diagrams in 
the decay B —>■ Xs'y containing heavy quarks and heavy gauge bosons in the loop. We 
will discuss the result for these diagrams, shown in Fig. [H along with all other new 
contributions in the next section. Concerning the notation, here and in the following we 
denote by V the heavy charged gauge bosons. 

In the Higgs sector the symmetry breaking structure leads to a number of charged 
and neutral Higgs fields, which in principle also transmit flavor changing interactions. 
Considering the heavy top quark coupling only, the quark sector can be identified with 
a modified 2 Higgs doublet model (2HDM) type II, in which only the third generation 
diagrams are taken into accountH. This corresponds to [8], where the Higgs sector has 
been approximated in this way. Note that these kind of terms do not appear in FCNC 
processes mediated by a tree level Z', since there the top Yukawa does not appear. 
However, for B — > Xs'y they should be added. 



3 Present Situation of the Decay B X57 

We begin with several general remarks concerning the decay B Xs'y, while the in- 
terested reader may consult [32] for a very recent and more elaborate discussion of the 
theoretical and experimental background. In general, theoretical interest in this decay 
stems from the following features: 

• Being an inclusive process, it can be calculated much more reliably than exclusive 
processes generally can be. Using the heavy quark expansion, the calculation can 

^Note also that due to the exchange of the top and bottom quark in the third quark triplet, the 
couplings are accordingly exchanged. 
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be performed at the partonic level, while additional corrections are suppressed by 
terms of (9(l/mb). 



• The decay is by now very precisely measured. In addition to the increasingly 
precise SM calculation, it therefore offers a very good test of the SM or a nice tool 
to search for physics beyond it. 

Theoretically the decay B — > Xg'y has now been calculated completely up to NLL preci- 
sion, and recently a first estimate of the value at NNLL precision has been obtained [33], 
after a great effort of many groups to calculate the different contributions required [34] . 
It reads: 

Br{B ^ X,7) Innll = (3.15 ± 0.23) x 10"^ (8) 

This effort of the NNLL calculation had become necessary, due to a large uncertainty 
of the NLL calculation resulting from uncertainties in the charm mass renormalization 
scheme and scale. 

On the other hand, the current experimental average is combined from the measure- 
ments of BaBar, Belle, CLEO and others, using inclusive and exclusive methods. Using 
a photon energy cut of > 1.6, HFAG gives an average of [42] 

BriB ^ X,7) Iexp = (3.55 ± 0.241°;?^ ± 0.03) x 10"^ , (9) 

where the errors include also those from the extrapolation to lower energies. Therefore, 
the experimental value is somewhat higher than the theoretical one, but still quite well 
compatible. 

Let us also note here that the decay B X^j has been analyzed in various models 
beyond the SM, among these the 2HDM type II [35] , the general as well as the minimal 
flavor violating (MFV) minimal supersymmetric SM (MSSM) [36], models with one [37] 
or two [38] universal extra dimensions (UED) and the littlest Higgs model with [39] and 
without [40] T parity. Some of these analyses have been performed at the NLO level, 
but most of them are done at LO order only. A comparison and compilation of the main 
results can be found in [32]. 

3.1 New Contributions to B ^ X,,-/ in the Minimal 331 Model 

In this section, we list all the different contributions to the decay B —>■ Xs'~f that modify 
the prediction of the 331 model with respect to the SM, and comment briefly on the 
cancellation of divergences. All expressions are explicitely given in [8]. Concentrating 
on the operators that are relevant for B Xs'~f, we have also confirmed the calculation 
of [8]. 

There are three different possibilities for the additional particles to show up in the 
decay B X^j. First of all, there are penguin diagrams involving the new charged 
gauge bosons, shown in the first two diagrams of Fig. [TJ For an arbitrary charged gauge 
boson and quark, the general structure, after calculating all penguin and the self energy 
diagrams required and summing over all of these, is a well known generalization of the 
Inami Lim Functions. In this calculation, a divergent term remains, which additionally 
violates gauge invariance. Within the SM, this term is canceled by the GIM mechanism 
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when summing over all quark flavors. In the 331 model, this singularity is not removed 
entirely even then, due to the different charge assignment of the quarks. In this case, the 



cancellation is achieved [8,9] by adding the Z' — 7 mixing diagram shown in Fig. 1(c) 
Therefore, the very same mechanism that generated the FCNCs at tree level can poten- 
tially also enhance them at loop level, in particular, if the GIM mechanism is as effective 
in the SM as it is in i? ^ ^s7- 



Next, there are also the Z' penguin diagrams shown in Fig. 1(d) To simplify the 



expressions, we assume that all quark masses vanish (in this case, there are down type 
quarks also in the loop). The unitarity of the mixing matrix V then allows to sum up 
the terms from the c?, c, and h quarks into a simple form without any leftover divergence. 

Finally, there are contributions from the Higgs sector of the model. Calculating these 
explicitely results in extremely complicated expressions, due to the elaborate structure 
of the Higgs mass terms. Resorting to the simplifications mentioned above, one can 
describe the Higgs sector by the 2HDM, where the corresponding expressions are well 
known [35]. 



3.2 Initial Conditions and Renormalization Group Analysis 

The standard procedure of calculating any weak decay is the framework of a weak effec- 
tive Hamiltonian. Here, the separation of scales is achieved by integrating out all heavy 
degrees of freedom, and casting the Hamiltonian into an effective form: 

Heff = -2^Gp C^{ii)0,{ii) . (10) 

i 

The sum extends over all operators that can appear at a given scale. Within the SM, the 
masses of the heavy particles are all of the order of the weak scale, at which initial con- 
ditions are then calculated. Next, the QCD corrections arising from the renormalization 
group running further enhance the branching fraction oi B ^ by about a factor 
of 3, due to the effectiveness of the GIM suppression in this case. Within the the 331 
model, the QCD corrections are not expected to be as important [8], but we take them 
into account for completeness. We use the leading order formulae for the RG running, 
with the anomalous dimensions as given in [8] , where they have been extended from the 
SM to include the additional operators present in the 331 case. Turning to the matching 
conditions in the 331 model, there are evidently now additional scales in the problem. 
In principle, there is the scale of the Z' boson mass, as well as that of the masses of the 
charged gauge bosons. The QCD running between these two scales does not modify the 
values of the coefficients by much, and we therefore integrate out all heavy particles at 
the lower scale My. These are, in addition to the heavy gauge bosons, also the heavy 
quarks, while the Higgs sector is added at lower energies. 

These initial conditions are run down to the weak scale, where several changes occur: 
First, the top quark is integrated out, and along with it the operators in which the top 
quark appears. These are replaced by the standard operators (5i/2 involving the charm 
quark. Next, we should here add the SM matching conditions as well as those from the 
2HDM. 
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Finally, from these coefficients, the branching ratio of -B — ^ X^'^ is calculated by 

in the notation of [41], where C = \Vub/Vcb\'^r{b -> ceVe)/T{b ueVe) = 0.58±0.016 and 
N{Eo) = 0.0031 are the nonperturbative corrections (we take Br{b — ceVe) = 0.106). 
Using the SM expression for Cf given above leads to the LO value of the branching 
fraction, which, numerically, does not agree with the most recent theoretical value given 
in [33]. To accommodate for the corresponding shift, we directly set the SM part of Cj 
to the NNLO value and go through the entire RGE procedure only for the new physics 
contributions. 



4 Numerical Analysis 

4.1 Preliminaries 

We will now analyze numerically the new contributions to the decay B — > and 
investigate whether additional bounds on the minimal 331 model can be obtained from 
this mode and which of the parameters appearing have the strongest influence. 

As constraints we will use the existing data from K and B meson mixing, such as 
AMd/s, sin 2/3, AM^ and Ek whose theoretical expressions in the 331 model are given 
in [6]. As numerical input we will take the tree level values of \ Vus\, \Vub\, given in [6], 
and 7 = (82 ± 20)°. Further, we follow the analysis of [39] and set all non-perturbative 
parameters to their central values and allow AMk, Sk, AMd, and AMg to differ from 
their experimental values by ±50%, ±40%, ±40% and ±40% respectively. In the case of 
AMg/AMd we will choose ±20% as the error on the relevant parameter is smaller than 
in the case of AM^ and AMg separately. Similarly, we will absorb the SM uncertainties 
of 5 — > Xs'~f, stemming mainly from the remaining scale uncertainty in the charm mass 
as well as the CKM factors, into the experimental uncertainty by increasing it to ±15% 
instead of the ±8% given above. 

We then perform a scan over the parameters of the 331 model considering the fol- 
lowing ranges: 

mr, ms = 250 - 1000 GeV, Mh+ = 300 - 2000 GeV, Mz' = 1000 - 5000 GeV. 

The three angles and the three phases of the new mixing matrix V are kept arbitrary. 
Further, the expressions for the 2HDM depend on tan/5, which we mainly restrict to 
values of tan /5 > 1 for reasons that will become evident during our analysis. We observe 
that, in a parameterization of the V matrix that keeps the real and imaginary parts of 
the relevant combinations of V*jVim as independent, the only bound we need to consider 
will be the one coming from AMg. 

4.2 Constraints from B Xgj Compared to Other Constraints 

Let us ffist elaborate on the influence of the chosen values for tan /5. Looking at Fig. [21 we 
observe that the value obtained for B Xg'j is practically independent of its numerical 
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Figure 2: The dependence of Br{B Xg^) on tan/3. For values of tan/3 below 0.5 significant 
enhancements are possible. The grey band indicates Br{B — > Xs7)|exp, with errors inflated as described 
in the text. The upper band shows the values for Af^+ = 300 GeV, while the lower one shows those 
for Mh+ = 5000 GeV 



value as long as tan/3 > 2. On the other hand, large values of the branching fraction 
can be obtained for smaller values of tan/3. This effect of the 2HDM is well known [35], 
and has most recently been investigated and numerically updated in [33]. We refer the 
reader to the detailed discussion oi B ^ within the 2HDM given there, and in 
the following fix tan/5 = 2, in order to show more clearly the additional effects of the 
331 model. In this context, we would like to point out that, in the pure 2HDM, very 
small values of tan/5 are excluded by other observables, such as electroweak precision 
tests. While the 331 model does no longer resemble the 2HDM when gauge bosons are 
included, we use only the larger values for tan (3. 

Next, we show in Fig. [3] the dependence of the branching fraction on the Z' mass, 
where we separate the dependence on the charged Higgs mass by showing only the values 
obtained for the upper and lower bound on Mh+ , respectively. In any case, the upper line 
corresponds to the lower value of Mh+, meaning that the branching fraction increases 
with decreasing Higgs mass. The width of the remaining bands corresponds then to the 
allowed range of V3*2V33 and the heavy quark masses. Clearly, a variation in the charged 
Higgs mass has a much stronger influence on the value of the branching fraction, which 
leads us to conclude that, as a whole, the 2HDM contributions vastly dominate over those 
from the gauge bosons, and that, therefore, the signature of the minimal 331 model with 
respect to B ^ Xsj is basically identical. On the other hand, the bounds on the Higgs 
sector from the type II 2HDM can immediately be applied to the 331 model. Also, we 
can conclude that the decay B Xj,7 would be well suited to explore the Higgs sector 
of the 331 model without additional pollution from gauge bosons, if the minimal 331 
model should be established through other channels. These results represent the main 
conclusions of our analysis. 

In this context, we also show the corresponding dependence of the branching fraction 
on the charged Higgs mass Mh+ in Fig. HJ which in principle allows us to read off a lower 
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Figure 3: Br{B Xg^) versus Mz' ■ The grey band shows the experimental range with a 15 % error. 
The upper band corresponds to a Higgs mass of = 300 GeV and the lower one to Mu+ — 2000 GeV. 

bound for Mu+ ^ 400 GeV. One should now take into account the NLO corrections 
within the SM and the 2HDM, which are known to be significant, in order to improve 
quantitatively on this bound. However, since B Xs'j has been studied elaborately, 
both in the SM and the 2HDM, in [33], including every known contributions to this 
decay, we consider it beyond the scope of our paper to repeat this analysis and simply 
quote the lower bound of the Higgs mass as 



Interestingly, while the 331 model can, in principle, either enhance or suppress the 
branching fraction of i? — Xs'j, depending on the sign of the mixing matrix elements, 
the model as a whole predicts an enhancement of the branching fraction, due to this 
strong dominance of the 2HDM model. This is, of course, rather fortunate considering 
the present experimental result. 

We will close this subsection and our numerical analysis with some very brief remarks 
concerning the possible infiuence on the rare decays and the other observables discussed 
in [6] . However, after the findings of the last paragraph, one expects this infiuence to be 
completely negligible. In any case, the most interesting observable is the Bg — B^ mixing 
phase, which can be large. This statement is not altered by our analysis of i? ^ Xg^. 
More interestingly, a first measurement of this quantity has recently appeared [52] , which 
now offers some hint as to its size. As of now, this measurement is not precise enough 
to warrant a more detailed discussion, but it will be extremely interesting to follow the 
progress of the corresponding experiments. 

All other statements of [6] remain unaffected, in particular, one can still obtain 
sizeable effects in the rare K decays, while large effects in Ba/s — > fi~^fi~ seem excluded. 
To show this graphically, we show in Fig.[5]the correlation between the branching fraction 
of -B — >• Xs'j with the ratio Br{Bs fi~^fi~)\33i/Br{Bs — > fJ'~^fi~)\sM, where we observe 
that the grey band, representing the experimental range of Br{B — >■ -^^7), still allows 
for the entire range of Br{Bs — > 



Mh+ > 295 GeV 
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Figure 4: Br{B ^s7) versus Mu+. The grey band shows the experimental range with a 15 % error. 
From the plot one finds a lower bound on the eharged Higgs mass of roughly 400 GeV. 
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Figure 5: Correlation oi B Xs"( with the deviation of Br{Bs) = Br{Bs ii^ji^) from the SM. 
The experimental range of the former (given by grey band) does not further constrain the range of the 
latter. 
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5 Conclusions 



Using new data for both the decay B Xg'^ as well as for the other existing constraints 
available, and in view of recent theoretical progress, including a first NNLO estimate, we 
have reinvestigated the implications of the minimal 331 model for the decay B X^^. 
In contrast to several other FCNC observables, that are affected at tree level by Z' con- 
tributions within this model, we are dealing here with a purely loop induced process, 
which, due to the breaking of the GIM mechanism, may still receive significant con- 
tributions. We have in general retained the general feature of the more than ten year 
old analysis performed in [8], which particularly concerns the performed simplifications 
within the Higgs sector. It is described here in terms of an effective 2HDM, so that the 
new contributions to the branching fraction can be discussed in terms of two different 
parts: one originating from the extended Higgs sector and the other one from the ad- 
ditional gauge bosons and quarks. The latter are governed by the new mixing matrix 
V , that appears as a set of new parameters in the model, and is constrained by the 
existing bounds from FCNC processes. On the other hand, the 2HDM contribution is 
only constrained by the usual parameter constraints on the Higgs mass. The main new 
result of our analysis is the finding that the gauge contributions are now constrained so 
strongly that effectively one is left with the 2HDM ones. Thus, it is possible to obtain a 
lower bound on the charged Higgs mass, as obtained from a recent, more sophisticated 
analysis of the 2HDM. We therefore conclude that the decay B —>■ X^'j is a very useful 
tool to probe the Higgs sector of the minimal 331 model and is therefore complementary 
to other FCNC observables. This shows once again the power of the B and K meson 
systems in obtaining information about models beyond the SM, when they are used 
together and combined. 
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